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The dispersion relation and damping rate of longitudinal color plasmons in finite strange quark 
matter (strangelets) are evaluated in the limits of weak coupling, low temperature, and long wave- 
length. The property of the QCD vacuum surrounding a strangelet makes the frequency of the 
plasmons nearly the same as the color plasma frequency of bulk matter. The plasmons are damped 
by their coupling with individual excitations of particle-hole pairs of quarks, of which the energy 
levels are discretized by the boundary. For strangelets of macroscopic size, the lifetime of the plas- 
mons is found to be proportional to the size, as in the case of the usual plasma oscillations in metal 
nanoparticles. 



PACS numbers: 12.38.Mh, 25.75.-q 



I. INTRODUCTION 



OS . At extremely high temperature and/or baryon density, quark-gluon plasmas Q are considered to be energetically 
favorable as compared with hadronic matter. Such extreme conditions prevent us from confirming the presence of 
quark-gluon plasmas, although it is expected in interiors of neutron stars and in ultrarelativistic heavy- ion collisions. 
The intriguing possibility that strange quark matter, composed of u, d, and s quarks, might be the ground state of 
| the strong interaction was pointed out by Witten Q). If this is correct, finite strange quark matter, which is usually 
referred to as strangelets, would be a more stable self-bound system than a 56 Fe nucleus. Then, disruption of neutron 
£"--» stars and highly energetic heavy-ion collisions could leave behind many strangelets. This might allow us to detect 

quark-gluon plasmas in the form of strangelets under various experimental situations ||-[5| ■ 
On ■ The ground-state properties of strangelets were described by Farhi and Jaffe || within the MIT bag model; finite- 
size effects such as surface, Coulomb, and shell effects were taken into account. In ultrarelativistic heavy-ion collisions, 
however, the possibly formed strangelets should be accompanied by thermally excited states. Among various exci- 
tations, we take note of spin- and flavor-symmetric, longitudinal color plasmons that represent collective coun- 
teroscillations of different color degrees of freedom. In a strangelet, these plasmons are expected not only to exhibit 
a dispersion relation dependent on the color dielectric property of the medium, but also to undergo a size-dependent 
damping as in the case of an optically excited metal cluster or nanoparticle. When the nanoparticle is small enough 
for the electrons to be nearly collisionless and for the electric dipole radiation to be of little significance but is too 
. i-h , large for the shell structure to take effect, the surface plasmon excitations are damped in a time proportional to the 
size • Their energy is resonantly absorbed by individual excitations of particle- hole pairs of the electrons having 
discrete energy levels due to the presence of the boundary. Notice that the quark-gluon plasma in a strangelet is fully 
relativistic in contrast to the electron gas in a nanoparticle. Consequently, new features should arise from spin-orbit 
interactions and from excitations involving antiquarks, both of which may affect the lifetimes of the color plasmons. 

Longitudinal color plasma oscillations in relativistic and degenerate bulk matter of u, d, and s quarks were inves- 
tigated ]t],H by calculating the current-current correlation function within the random-phase approximation (RPA) 
which is valid in the weak-coupling limit. The dispersion law and damping for these oscillations are essentially the 
same as those examined by Jancovici jllj for electromagnetic plasma oscillations in a relativistic and degenerate elec- 
tron gas. In either case, the mode properties are characterized by the corresponding plasma frequency and Landau 
damping as in a nonrelativistic plasma. In this paper, we estimate the long-wavelength, low-temperature properties 
of longitudinal color oscillation modes in a strangelet of macroscopic size, surrounded by the QCD vacuum, by in- 
corporating the discrete energy levels of quarks into the current-current correlation function within the RPA. The 
vanishing color dielectric constant of the vacuum leads to the absence of the surface color charge; if present, it would 
modify the effective color electric field inside the strangelet according to a Mie-type formula. We thus find that the 
frequency of the longitudinal color plasmons in the strangelet is almost the same as that in bulk quark matter. The 
lifetime of the plasmons is determined by the dipole transitions from a quark state in the Fermi sea to an unoccupied 
one. It is shown that the quark particle-hole pairs involved damp the plasmons at a rate inversely proportional to the 
strangelet size. Beyond the RPA, however, the plasmons are further damped by collisions between quarks. 

In Sec. II, the longitudinal polarization function is calculated up to one-loop order, allowing for the discrete quark 
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states inside a strangelet. The resulting dispersion relation and damping rate of the longitudinal color plasma oscil- 
lation are described in Sees. Ill and IV, respectively. Concluding remarks are given in Sec. V. 



II. POLARIZATION 



We consider a spherical strangelet of radius R and baryon number A, composed of a nearly noninteracting and 
uniform gas of u, d, and s quarks and embedded in the QCD vacuum. The chemical potential /i.; of quarks of flavor i, 
including the rest mass m,, is assumed to be independent of i. Such independence ensures flavor equilibrium among 
the quarks, since normally leptons have a vanishingly small chemical potential in the strangelet ||. We work in units 
in which Ti = c = ks = 1- We set /i,; = /i = 300 MeV, m u — md — 0, and m s = 0-300 MeV. The temperature T of the 
system is taken to be in the range < T -C [i. The zero-temperature model for quark matter will thus be adopted 
below. 

Let us now examine the properties of longitudinal color oscillation modes in a strangelet by calculating the cor- 
responding proper polarization IIoo in the RPA (one-loop order). At T = 0, IIoo can be expressed in terms of the 
correlation function of the time component of the quark current (a = 1, . . . , 8) as 

Iloo(q) = -y // x j/y / d(x°-y )e^'^{^)\Mx),J°(y)])e(x -y°) , (1) 

with 

JZ{x) = ASb&i(aOYVw(aO > 

where q = (w,q) is the four-momentum transfer, V = 4irR 3 /3 is the volume of the strangelet, g is the renormalized 
color coupling constant, tp ai is the quark field of color a and flavor i, \™ b are the SU(3) generators, and (■ • •) denotes 
the expectation value in the ground state. The right-hand side of Eq. (1) is not summed over a and is invariant under 
change in a. In the limit V — > oo, this form of Hoo(q) reduces to the usual form written for bulk uniform matter. 
Hereafter, we assume that the momentum transfer |q| is far smaller than R . 

The essential part of the present study is to incorporate the discrete quark levels into the one-loop order calculations 
of n o- For such calculations, we first write the quark field in the form of a mode expansion 

<M*°,X) - Y^teKai^lM)^* + 4-^( X ) e ^*°] > ( 2 ) 



where ^^ ai and c^ Ka i are the spinor and the annihilation operator for the quark eigenstate of energy LUki = \/fc 2 + mf 
characterized by a set of the quantum numbers £ = {k,j,m}, with wave number k, angular momentum j, and its z 
component m, and by the integer k = ±(j + 1/2) corresponding to the parity (— 1) 3±1 / 2 ; ~&£ Kai and d^ Ka i are for the 
antiquark eigenstate of energy u>ki- The spinors ^^ ai can be obtained from the Dirac equation 

(<7 • V ± 7°w w - rn^fy = for r < R (3) 

and the boundary condition 

- *7-?*£L = *gL for r = R . ( 4 ) 

where r and f are the length and the unit vector of the position vector x, of which the origin is set to be the center of 
the strangelet. Condition (4) ensures the absence of the outgoing probability flux from the boundary. The analytic 
solutions to Eq. (3) under condition (4) were written |12|,[l3[ in the context of the MIT bag model; we adopt the 
solutions satisfying the normalization condition Jy dV (^^^^^^ — 1. Such solutions read 



with 



9ltl(r) = A<£>j K (kr) + B^n K {kr) , (6) 
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ftliT) = — — [A ( ^i(kr) + Bg^W] . (7) 

±Ulki + TTli 



Here % a is the color eigenstate, j n and n n are the nth-order spherical Bessel functions, 

1-1/2, 



=J + 1 /2,oA - 



y/j-m + lY£-;<'(Sl) 



„ jm («) = —— _^^;i/2_ ) for« = j + l/2 (8) 



and 



(9) 



are the two-component spherical spinors, and 



kR 

; {±2u kl R{K±uj kl R) + m l R]- 1/2 , K = j + l/2, 



and 



4± } = ^ (fci?)|i?3/ (io) 

lo, k = -Cj + 1/2), 

f , k = j + 1/2 , 

<1= jfci? 1/9 ( n ) 

1 u_ ^7fcWR372 [ ±2 "fc«fl(« ± "MR) + m*R3~ 1/2 , K=-(i + l/2) , 

are the normalization factors. The energy eigenvalue E = ±uj k i may be determined from the resultant quantization 
condition ]l4| 

3l=j±i/2(kR) = T E + ji=3TiM kR ) for K = ±(j + i) . (12) 

In the case of interest here, in which the baryon number A is macroscopic, the solutions to Eq. (12) approximately 
read k = Niir/R with the principal quantum number Nj = 1,2,3,... . Here, the asymptotic form of jiikR) ~ 
(1/kR) cos[fci?— (I + 1)tt/2] for large kR has been used, and only quantities of leading order in kR have been retained. 
The JVjth solution has a degree of degeneracy, grjVj = 12(Ni + l)(2Ni + l), stemming from the color, angular momentum 
j = 1/2, 3/2, . . . , 2kR/n — 3/2, and parity (— 1) J±1 / 2 [M. We find from this that the Fermi wave number of i 
quarks is kp ri rs (7r 4 n i /6) 1 / 3 with number density ri,, which is in turn related to [ii as fa = Lu kFii . 

By building the resulting quark field t/j al into Eq. (1) and by noting (c| Kai c 5KM ) = 6(fj,i -ojm) and {d\ Kai d^ Ra i) = 0, 
calculations of n o up to lowest order in q can be performed straightforwardly. The result, being of order q 2 and 

expressed in terms of dipole transitions, can be formally divided as Uoo(q) = Il^ at '{q) + U^\q). IIqq 121 ''' is the 
matter contribution given by 



with 



x ( —■ * 7^ L — ) M tl^M) + ~ ^i) 

\UJ + IT] + LU kl i - UJk 2 i W + IT] - UJkii + 0Jk 2 i J ? S 

x [ _ - \ M {m) fall (13) 

\uj + n] + uj kll + uj k2l lu + n] - Lu kll - u k2 J *i&'W H 'J » v ; 

Ml*L 1K2 = ~ //* //y[q-(x-y)] 2 *wj^ 

q 2 ,/, . . . . . , V|2 \{h-mi + l)(ji + mi + l) 
= -777\{ki,Ji,n 1 ,i,+\r\k2,j2,K 2 ,i,+)\ 



(Ji + I) 2 

„x x at \ , C7l +mi) . , 

mi ,rri2 

Jl 

o 



1 r r /i / \ 

-2/- T 1 \2 Jl ,h"mi,m2"\~ K1K2) 
Jl Ul X J 



(14) 
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M tL^ = ~w J dx I d y^< x -y)] 2 *SL i ( x )*Ua 2i ( x )^ 2 i(y)*^a li (y) 



(J i - mi + l)Q'i + mi + 1) 
(71 + I) 2 



— |(fci, ji,Ki,i,+|r|fc 2 , j 2 ,K 2 ,i,-)\ 2 

X X Of \ , (il - m l)0'l + m r r /,/ \ 

X0 il+l,i2 () mi,"i2 fc 'l K l K 2j H 3 — 1 , j2 "mi ,m 2 "v'^l'*2 J 



J1U1 + 1 ) 



and IIqq P '' is the vacuum contribution given by 



n& p) (,)--2(f) 2 E E 



* £l?2KlK2 



1 



1 



lu + in + uj kli + uj k2l ui + in - uj kl i - u k2l 



M 



(99) 



ClC 2 KlK2 



(q)- 



(15) 



(16) 



Here 77 is a positive infinitesimal, and (fci, ii, Ki, i, +|r|fc2,j2, «2, ±), in which the sign + (— ) specifies the quark 
(antiquark) state, denote the radial components of the matrix elements. These components vanish when k\ — fc 2 , or 
else they are given by 



|<fci,ii, /si = ±(ii + 1/2), i, +\r\k 2 ,j 2 = j 1 + 1, k 2 = ±(ii + 3/2), i, +) 



2fcifc 2 



(fc 2 - klfR 



[{ljJ kl j + LO k2 j) 2 R 2 + («i + AC 2 )(^fc 1 j + ^fc 2 j)fi] 2 

[2wfe 1 ji?(Ki + uj kll R) + m l R][2uj k2i R{K 2 + uj k2l R) + m l R] 



(17) 



fci.ii, ki = ±(ii + 1/2), i, +\r\k 2 ,j 2 = ii - 1,K2 = ±(ii - 1/2), i, +)| : 



2fcifc 2 



(fc 2 - fc 2 ) 2 i? 



[(q; fcl i + uj k2 i) 2 R 2 + («i + K 2 )(u fclt + ^ fc2t ).B] 2 
[2w fel ji?(Ki + Wfc^i?) + mji?][2w fe2i i?(K 2 + oj k2l R) + m l R] 



(18) 



|(fci,ii, ki = ±(ii + 1/2), *, +|r|fc 2 , j 2 =ii,K 2 = T(ii + l/2),i,+)|' 

i 2 



2mi 



u kl i — ljJk 2 i 



+ («1 - Ki) 



/ 2hk 2 R \ 2 _ 

V k\ - k\ J [2Lu kli R(ni + Lu kli R) + rriiR] [2u> k2i R(n 2 + u k2i R) + m,i?] 



(19) 



Kfci.ii, /si = ±(ii + l/2),i,+|r|fc 2 ,i 2 =ii + 1,K2 = ±(ii+3/2),z,- 



2fcifc 2 



(M - fc 2 2 ) 2 i? 



[(^fcii ~ u k2l ) 2 R 2 + (ki + K 2 ){u kl i - ^fc 2 »)-R] 2 
[2w feli i?(Ki + u> kli R) + miR][2uj k2l R(-K 2 + uj k2l R) + rriiR] 



(20) 



|(fci,ii, /si = ±(ii + l/2),i,+\r\k 2 ,h=ji - 1,K2 = ±(ii - 1/2), i,-) I 5 



2fcifc 2 



(ki - kl) 2 R 



[((J kl i - Uk 2 i) 2 R 2 + («1 + K 2 )(w fel4 - Wfe 2J )i?] 2 



[2w feli i?(«;i + w feli .R) + miR]{2u; k2l R(-K 2 + u k2l R) + niiR] 

|(fci,ii, Ki =±(ii + l/2),i,+|r|fc 2 ,i2=ii,K2 = T(ii + l/2),i,-)| 2 
_ / 2kik 2 R \ ' 





(u) kli +uj k2i )R + 


2mi 


+ («i - K2) 


2 








[2u) kli R(ni + uj kli R) + m 


i?][2^ fc2J i?( 


-K 2 + Wfc 2i i?) + niiR] 



(21) 



(22) 



It is convenient to further divide the matter part n o ™ at ^ as n o ™ at ^ = ilg^ + + , where 11^ , I1 q ^ , and 
n o ^^ are composed of the matrix elements written by Eqs. (17) and (18), Eq. (19), and Eqs. (20)-(22), respectively. 
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We thus find that H oc j (n o g ) arises from excitations of quark particle-hole pairs, accompanied by the change in the 

angular momentum of the scattered quark Aj = j% — ji = ±1 (Aj = 0); IIqq 9 '' (1T q P ' > ) arises from excitations of 
quark-antiquark pairs via scattering of the quark in the filled Fermi (Dirac) sea, accompanied by Aj — 0,±1, onto 

the vacant Dirac (Fermi) sea. The quantities Hqq , n Q 9 \ and II q P ^ are associated with effects of relativity: the 
spin-orbit interactions, fluctuations of the Fermi sea involving quark-antiquark pairs, and the second-order vacuum 
polarization, respectively. 



III. DISPERSION RELATION 



The dispersion relation of the longitudinal color plasma oscillation has been obtained from the one-loop polarization 
IIoo given by Eqs. (13) and (16). By noting the complete analogy with the case of a metal nanoparticle in a dielectric 
medium jicj , such a relation can be written as 



1 Relloo (g) 

7T £ t 

3 q 2 



, 



(23) 



where eo is the color dielectric constant of the strangelet, having the color carrier part subtracted out, and e vac is 
that of the surrounding vacuum. We normally set eo = 1 and e vac = 0. Up to leading order in kp.iR, the summation 

over £i, £2, and ^2 m Eqs. (13) and (16) may be performed by, if available, replacing Y] . — > Jq 2 *^" dj and 



J2k ~ * (R/ 77 ) Jo F * ^he mam contribution to Relloo is ReII 0C / = q 2 w 2 /3u 2 , where u p — 9 2n i/^^i is the 

color plasma frequency of bulk quark matter Q^]. |ReII 



l/kp ti R (a s ) with respect to |RcIT Q' ) |/q 2 (unity), where a s — g 2 /4w is the QCD fine structure constant. We thus 
obtain the mode frequency u> ~ lo p . Its agreement with the bulk value stems from absence of the surface color charge 
due to e vac = 0. The values of u> p calculated as a function of rj s = m s j [i for a s = 0.4 have been exhibited in Fig. 1. 
We have found that uj p depends on m s only weakly. 



(so) 

00 



l/q 2 (|Ren 00 



ReliQQ^I/q 2 ) is suppressed by one power 



IV. DAMPING RATE 



Within the framework of the RPA, the damping of the longitudinal color oscillation mode is determined by its 
resonant energy transfer into individual excitations of particle-hole pairs of the quarks whose eigenenergies are made 
discrete by the boundary. The corresponding damping rate may be derived from the polarization IIoo obtained above 

as a 



Trpa = -■^ImIIoo(<7)| [ , 



(24) 



Here Trpa <C w p is assumed, and the imaginary part of IIoo (?) is expressed up to leading order in kp^R and q as 



Imlloo(g) 



T7 4 /U,R 



o 2 E 



h(u,T]i) 



(25) 



with u = u)//j,, rji = rtii/ n, and 



H^vd = - 



dz 



" J max(l— v,y)i) 



\/z 2 -rjjy/(z + v) 2 -rji 
(2z + v) 2 



/9 9 , -k 2 v(%z 3 + 8vz 2 + Av 2 z + V 3 ) , 
64z(2z + v) 



2 ,2 



n 2 v(&z 3 + \Qvz 2 + Ylv 2 z + 3is 3 ) 
M(z + v)(2z + v) 



ln 



1 - 



* - vi 



z~ - m 



7T J (Z + Vf 



Z2z{z + v) 



1 + ln 



l + -\l z -m^R 



(26) 
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In Fig. 2, we have shown Trpa/wj, as a function of R for various values of r) s . For lo < fi + rrii, which is satisfied at 

u> = ujp, ImiTQQ P ' > cancels out ImHQQ 9 ''. On the other hand, the contributions of Inin^ and Inin^ ', the latter being 
suppressed by one power a s with respect to the former, lead to the damping of which the rate behaves as ~ c/R 
|l6| . Such contributions may be regarded as Landau damping, i.e., a decay of the collective mode by exciting a single 
real particle-hole pair. In this case, the surface of the strangelet, which acts as a boundary condition determining 
the discrete quark levels rather than as a scatterer, renders Trpa nonzero quantum mechanically and independent 
of a s in the weak-coupling limit. As can be observed from Fig. 2, Trpa has only a slight dependence on m s . In 
the nonrelativistic limit, only the part Indl^ remains up to 0(l/kp^R), still yielding Trpa oc 1 . We remark in 
passing that for a strangelet with radius R < 3 fm, which has a baryon number of less than ~ 30, higher-order terms 
in (fc^ii?) -1 as ignored here may well have consequence to the frequency and damping rate of the color plasmons. 

We next consider the damping of the longitudinal color plasmons due to collisions |l7[ between quarks in the weak- 
coupling regime. These collisions, by color exchange, disrupt the self-consistent fields responsible for the plasmons 
and add r co i to the RPA result Frpa- If one attempts to calculate r co i in perturbation theory, one meets an 
infrared divergence that arises from the color magnetic interactions mediated by transverse gluons. At least for finite 
temperature fl§|| , such a divergence cannot be removed even by dynamical screening fl9| of the magnetic interactions 
p0| . This indicates that calculations of r co i up to leading order in a s require nonperturbative features of the relativistic 
quark-gluon plasma which may possibly develop a magnetic mass m mas and thereby provide an effective cutoff [^l| . 
For a thermal, massless plasma at zero baryon density, where m mag ~ a s T, one thus obtains r co i ~ a s ln(l/a s )T from 
kinetic theory jl7| . This rate arises from the collision integral associated with one-gluon exchange in the bulk plasma; 
finite-size corrections of 0(l/kF,iR) to such a bulk result for r co i should be of higher order in a s than Trpa, which 
remains in the absence of the collision integral. At the nonzero strange quark mass, baryon density, and temperature 
of interest here, however, the values of m mag and hence of r co i are essentially unknown. Therefore, all we can mention 
is that such a collisional damping or the Landau damping due to individual excitations of quark particle-hole pairs 
dominates the total damping rate r tot = Tp>pA + r co i, according to whether the mean free path of the quarks ~ 
is small or large compared with R. Here, it is instructive to note that no damping of the plasmons due to gluon 
radiation occurs in the absence of the surface color charge. 



V. CONCLUDING REMARKS 



We have calculated the dispersion relation and damping rate of the longitudinal color plasma oscillations in a 
strangelet in the regimes of weak coupling, low temperature, and long wavelength. For kp^R 3> 1, we have found 
that the frequency of such oscillations is almost as large as that of bulk color plasmons and that the Landau damping 
yields a lifetime proportional to R. However, the present analysis is accompanied by several uncertainties. We have 
confined ourselves to the case of low |q| and T. With increasing |q|, multipole plasmon excitations should develop. 
At finite temperature, furthermore, gluonic excitations in the intermediate state contribute to the proper polarization 
n o; within the one- loop approximation, such a contribution is simply added to the quark- loop contribution without 
interference ||. A study allowing for those excitations has yet to be done. Shell effects in strangelets also remain to 
be taken into account; such effects are predicted to manifest themselves remarkably for A < 100 j^,|2|. It is expected 
from the analogy with the case of collective excitations in normal nuclei that whether strangelets are in the vicinity 
of closed shells or not plays a role in the properties of the color plasma oscillations. Third, it is to be noted that 
the present analysis is restricted to the weak-coupling regime (a s < 1). In the strongly interacting nonperturbative 
regime (a s ~ 1), more relevant to a real system, the damping rate of the color plasmons could be greatly modified. 
These nonperturbative effects should have primary significance to the collisional damping. 

To conclude, we briefly address what type of traces a strangelet, if formed in heavy-ion collision experiments |23j ] 
performed at a beam energy of order 200 GeV per nucleon on a fixed target, could leave via excitations of the 
longitudinal color plasmons. Immediately after central nucleus-nucleus collisions at such energies, it is expected that 
a quark-gluon plasma would be produced in the resulting matter of sufficiently high energy density. Then, the quark- 
gluon plasma might develop into strangelets via pre- freeze-out evaporation of pions, K + , and K° which carry away 
entropy and antistrangeness from the system. This strangeness distillation mechanism, suggested by Greiner et al. 
p4j, would work as long as strangelets are stable against strong decay. Note that the freezeout temperature ~ 100- 



200 MeV 25 1 is comparable with the frequency ui p for realistic values of a s in the range ~ 0.4-1. Longitudinal color 
plasma oscillations would thus be thermally excited in strangelets until the temperature fully drops. On a time scale 
of order T^, the excited color plasmons would decay into quark particle-hole pairs, leading to emission of photons and 
lepton pairs via electromagnetic interactions. These electromagnetic signals, having energy of order or less than ~ u pi 
might be observed in the highly energetic heavy-ion collisions at CERN's SPS experiments ]23[ |. In order to support 
the picture mentioned above, however, estimates allowing for relaxation processes occurring in such a cooled and 
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decompressed system, for other mechanisms producing photons and lepton pairs, and for the uncertainties discussed 
in the preceding paragraph, are required. Elucidation of this problem is beyond the scope of this paper. 
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FIG. 1. Color plasma frequency lu p as a function of rj s , calculated for fi — 300 MeV and a s = 0.4. 
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FIG. 2. Damping rates Frpa divided by u) p as a function of the strangelet radius 7?, which have been evaluated for fi = 300 
MeV and a s = 0.4. The lines from bottom to top are the results obtained for rj s = 0, 2/3, 5/6, 1. 
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